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ABSTRACT 



A near-infrared excess is detected at the white dwarf PHL5038 in UKIDSS photometry, consistent with the presence of a cool, 
substellar companion. We have obtained H- and A"-grism spectra and images of PHL5038 using NIRI on Gemini North. The target is 
spatially and spectrally resolved into two components; an 8000 K DA white dwarf, and a likely L8 brown dwarf companion, separated 
by 0.94". The spectral type of the secondary was determined using standard spectral indices for late L and T dwarfs. The projected 
orbital separation of the binary is 55 AU, and so it becomes only the second known wide WD+dL binary to be found after GD 165 AB. 
This object could potentially be used as a benchmark for testing substellar evolutionary models at intermediate to older ages. 

Key words, white dwarf - brown dwarf - resolved binary 



1. Introduction 

The study of low mass companions to white dwarfs (WDs) 
probes high mass ratio binary formation and evolution, includ- 
ing the known deficit of brown dwarf (BD) companions to 
main-sequence stars (McCarthy & Zuckerman 2005; Grether & 
Lineweaver 2006) . There are currently no empirical constraints 
on substellar evolutionary models at intermediate to older ages 
(> lGyr), such as those expected for most WDs. Thus wide, de- 
tached substellar companions to WDs can provide benchmarks 
for testing BD cooling models at these ages (Pinfield et al. 
2006). Furthermore, the Earth-sized radii of WDs provide excel- 
lent contrast for the infrared detection of these cool, self lumi- 
nous companions (Probst 1983a, b; Zuckerman & Becklin 1987; 
Burleigh, Clarke & Hodgkin 2002; Farihi, Becklin & Zuckerman 

2008) . 

However, BD companions to WDs are rare (< 0.5%; Farihi, 
Becklin & Zuckerman 2005). So far only 3 L-type compan- 
ions to WDs have been discovered; GD165B (L4, Becklin & 
Zuckerman 1988; Kirkpatrick et al. 1999), GD1400B (L6-L7, 
Farihi & Christopher 2004; Dobbie et al. 2005; Burleigh et al. 

2009) and WD 0137-349B (L8, Maxted et al. 2006; Burleigh et 
al. 2006). Of these, GD165 is the only spatially resolved widely 
separated system. 

Here we report a near-infrared (NIR) excess at the WD PHL 
5038. Follow up observations reveal a spatially resolved substel- 
lar companion of late L-type; only the second such pair found 
(over 20 years later) after GD 165 AB. 



2. Target Identification 

PHL 5038 was first reported by Haro & Luyten (1962) in a pho- 
tographic search for faint blue stars near the South Galactic 
Pole with B — 18 mag, and a poorly constrained position of 
22 h 20.5 m - 00°43' (J2000). This star was recovered in the Sloan 
Digital Sky Survey as SDSS J 222030.68-004107.3 (Eisenstein 
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Fig. 1. Near-infrared acquisition images of PHL 5038 taken with 
Gemini + NIRI immediately prior to H- and A'-grism spec- 
troscopy. Frames are 6" across with north up and east left. The 
brown dwarf companion, PHL 5038B, is located 0.94" from A at 
position angle 293.2", as measured in the /f-band image shown 
above. 



et al. 2006) , a WD with r eff * 8037 K and logg = 8.28, yielding 
an accurate position for the WD. 

The WD was detected in the UKIDSS Large Area 
Survey (LAS), with DR4 photometry for PHL 5038 listed in 
Table 1. These values give an excess of 7 - 10cr when com- 
pared to the predicted values for the WD of (J,H,K) - 
(16.86, 16.72, 16.72)mag, calculated by folding our WD model 
(Section 4.3) through the appropriate UKIDSS filter response 
curves. 

The UKIDSS data flag (mergedClass) indicated the possibil- 
ity that PHL 5038 was an extended source at one or more band- 
passes. Follow-up observations were warranted to investigate the 
nature of the excess flux and its possible spatially extended na- 
ture . Thus, this object was included in a Gemini programme to 
clarify the nature of NIR excesses to WDs in the UKIDSS LAS. 
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Table 1. Near-Infrared photometry of the PHL 5038 binary system and nearby calibration star. 



Bandpass 


A 


PHL 5038 


PHL5038A 


PHL5038B 


2MASS J22203205-0040210 




(jim) 


(mag) 


(mag) 


(mag) 


(mag) 


Y [ 


1.03 


16.89 ±0.01 








./' 


1.25 


16.73 ± 0.02 








H [ 


1.63 


16.44 ± 0.03 








r 


2.20 


16.31 ±0.04 








H 1 


1.65 




16.82 ±0.07 


17.84 ± 0.08 


14.447 ± 0.053 


K 2 


2.19 




16.72 ±0.07 


17.18 ±0.08 


14.161 ±0.051 



'UKIDSS passbands as described in Hewett et al. (2006). 2 H- and ASband photometry were obtained on a single night. The photometric errors were 
calculated using a 5% error from the nearby 2MASS star used for flux calibration, a 2 - 3% photometric measurement error, and an assumed 5% 
internal error typical of near-infrared arrays. The final column displays the nearby 2MASS calibration star with its 2MASS catalogue photometry 
and errors. Colour transformations between the filters are of the order of 1% (Carpenter 2001) and so are neglected. 



3. Observations 

Near-infrared grism spectroscopy of PHL 5038 was obtained at 
Gemini North using NIRI (Hodap p et al. 2003) . During the first 
acquisition images, the science target was spatially resolved into 
two components separated by 0.94", and the slit was placed over 
both point sources for all subsequent programme observations. 
On 24 Aug 2008, a single 30 s //-band acquisition image was 
obtained prior to placement of the science targets on the slit. 
Spectroscopy was performed with the H grism and the 6 pixel 
slit (covering 1 .43 - 1 .96 fj.m at R « 500) with eight, single coadd 
frames of 225 s each at one of two offset positions along the 
slit, for a total exposure time of 30 min. The A0V telluric stan- 
dard HIP 102631 was observed with eight frames of Is and three 
coadds. Internal fiat field and arc lamp calibration frames were 
taken immediately following the science target spectroscopy. 

Three attempts were made to obtain Zf-band spectroscopy 
of PHL 5038. On 12 Aug 2008, the science target was imaged 
for 30 s at two offset positions, then placed on the slit, but the 
spectroscopic sequence was abandoned before usable data were 
taken. However, the acquisition image from this date could still 
be used for photometry. On 21 Aug 2008, the target was also 
imaged prior to being placed on the slit. Spectroscopy was per- 
formed with the K grism and the 6 pixel slit (covering 1 .90-2.49 
fim at R w 500), with six, single coadd frames of 225 s each 
at one of two offset positions along the slit, for a total integra- 
tion time of 22.5 min. The A0V telluric standard HIP 102631 
was observed with eight frames of 1 s and three coadds. Internal 
flat field and arc lamp calibration frames were taken immedi- 
ately following the science target spectroscopy. Unfortunately, 
the fainter of the two point sources was poorly aligned with the 
slit, resulting in a significant loss of spectroscopic signal for that 
component (see below). On 6 Nov 2008, an identical Zf-band 
imaging sequence preceded a successful placement of both point 
sources on the slit. However, due to poor seeing the system was 
not spatially well resolved, and the S/N was too low to merit an 
extraction for PHL5038B. 



4. Data and Analysis 

All NIRI frames were processed using Gemini IRAF niri and 
gnirs packages, including the creation of appropriate flat-field 
images, bad pixel masks, and wavelength calibrated spectral 
data. 



4.1. Photometry 

Since only a single //-band image was obtained, this frame 
was cleaned of bad pixels and flat fielded. For the Zf-band im- 
age pairs, the frames were cleaned of bad pixels, flat fielded, 
pair subtracted, shifted and recombined. Fully reduced H- and 
Zf-band images are shown in Figure 1; the WD is designated 
PHL5038A and the cool companion is PHL5038B. 

Aperture photometry and astrometry on these images was 
performed with the IRAF task apphot using aperture radii of 2, 
3, and 4 pixels (0.1 17" pixeL 1 ) and a sky annulus of 16-20 pix- 
els. For most frames, photometry obtained with the 4-pixel aper- 
ture radius was discarded owing to probable contamination from 
the other binary component. A single H- and Zf-band measure- 
ment was made for both PHL 5038 A and B. The science target 
photometry was calibrated using a nearby (52" distant) field star, 
2MASS J22203205-0040210, which fell onto the NIRI array in 
all imaging observations. Photometric data for PHL 5038 A and 
B, as well as the calibrator are given in Table 1 . 

4.2. Spectroscopy 

The spectroscopic frames were cleaned of bad pixels, flat 
fielded, pairwise subtracted, shifted, combined, and wavelength 
calibrated. A slightly more aggressive (relative to the default 
nscombine settings) outlier rejection algorithm was applied in 
the next to last step, which produced visibly superior spectral 
images for extraction. In the //-grism observations, the spectra 
of PHL 5038 A and B were well resolved spatially, and each was 
extracted using a 2-pixel spectral aperture. For the Zf-grism spec- 
tral extraction, we only used the data from 21 August 2008. In 
these observations the two science target spectra were spatially 
resolved, but the signal of component B was significantly fainter 
than expected due to poor slit alignment. The data are still us- 
able and a spectrum of each component was extracted using a 2- 
pixel spectral aperture, but care had to be taken with PHL 5038B, 
and the default aperture function had to be carefully corrected to 
avoid automatically picking up the signal of the WD. The telluric 
standard data was reduced and extracted identically. 

Flux calibration and telluric feature correction were per- 
formed for each order with Spextool (version 3.3; Cushing, 
Vacca & Rayner 2004; Vacca, Cushing & Rayner 2003), using 
the general xtellcor package. Although in principle this should 
produce a well calibrated spectrum, there are a couple of reasons 
why this proved difficult or impossible in this case: 1) the telluric 
standard and science observations differed by 0.2 airmasses, and 
2) slit losses; both ordinary and those caused by clear star-slit 
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Fig. 2. Modelled (black) and observed (grey) HK spectra of PHL5038 A plotted with photometry as measured in the acquisition 
images. Note the detection of Pa ff absorption. 




Fig. 3. Observed HK spectrum (grey) of PHL 5038 B plotted with photometry as measured in the acquisition images. An L8 spec- 
trum (2MASS 1632+1904) has been scaled to match the flux of the secondary and over-plotted (black) for comparison. 



misalignment. Therefore the spectral data were calibrated from 
the Table 1 photometry; some inter-order adjustment was neces- 
sary for both the WD and the companion. 

The spectral data were further processed to interpolate over 
data points where residual bad pixels or cosmic rays persisted 
(around several to one dozen), to crop the noise dominated 
edges, and to merge the data into a single, contiguous spec- 
trum for each binary component. The fully processed spectra for 
PHL 5038A and B are shown in Figures 2 and 3 respectively. 

4.3. Spectral Analysis 

An atmospheric fit of the SDSS spectrum of PHL5038A us- 
ing the models of Koester et al. (2005) yielded values of T eS = 
8000 + 100K and logg = 8.2 + 0.2, implying a mass of M W d = 
0.72 + 0. 15 M© and a distance of 64+ lOpc. The model spectrum 
is plotted over the observed spectrum in Figure 2 and is consis- 
tent with the detection of Pa„ absorption. No contamination from 
the secondary can be seen indicating we have successfully sepa- 
rated the components in our spectral analysis. 

To determine the spectral type of PHL 5038B, we measured 
the standard spectral indices for late L dwarfs and T dwarfs 
(Burgasser et al. 2006), and listed them in Table 2 with estab- 
lished L8, L9 and TO indices for comparison. The errors were 
determined by taking the standard deviation over the spectral in- 
dex denominators for the H- and TT-band regions. These values 
indicate the spectral type is most likely L8-L9, with the H2O-H 
index suggesting a spectral type of TO is less likely . 

The observed HK magnitudes (Table 1) are also consistent 
with a late L-type BD at a distance of 64 pc. Although the S/N 
of the spectrum is insufficient for the detection of any distinct 



absorption lines, the overall shape of the spectrum is also in 
best agreement with a spectral type of L8. The spectrum of an 
L8 field BD (2MASS1632 + 1904, McLean et al. 2003) has 
been scaled to 64 pc and plotted over the observed spectrum of 
PHL5038B in Figure 3 for direct comparison. 

5. Discussion 



Caballero, Burgasser & Klement (2008) list the spatial density 
of L8-L9 field dwarfs as 0.64 x 10 3 pc 3 . It follows that the 
number of expected L8-L9 field dwarfs within a cylinder of an- 
gular radius =1" and length = 20 pc centered on PHL 5038 A is 
» 3 x 10~ 9 . Therefore, PHL5038A and B are almost certainly a 
physical pair. 

The 0.94" angular separation of the binary, at the esti- 
mated distance of 64 pc equates to a projected orbital separation 
of 55 AU. When the main-sequence progenitor of PHL 5038 A 
evolved into a WD the orbital separation of PHL5038B ex- 
panded by a maximum factor of M M s/M WD (Jeans 1924; 
Zuckerman & Becklin 1987; Burleigh, Clarke & Hodgkin 2002). 
For M W d =0.72 M Q and M M s =3.2 M , this factor would have 
been « 4.4. Thus, the original projected semi-major axis was 
>13AU. 

The only other resolved L-type companion to a WD cur- 
rently known is GD 165B, which has a projected orbital sep- 
aration of 120 AU. Whether or not the secondary in this sys- 
tem is a true BD is still a matter of some debate (Zuckerman & 
Becklin 1992; Kirkpatrick, Henry & Liebert 1993; Kirkpatrick 
et al. 1999; Zuckerman et al. 2003). We estimate the temperature 
of PHL 5038B to be 1400- 1500K by comparison to L-type field 
dwarfs with measured effective temperatures (Vrba et al. 2004). 
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Table 2. Spectral indices as measured for PHL 5038 B, L8, L9 and TO fields dwarfs (Burgasser et al. 2006). 



Index 


PHL 5038 B 


2MASS 1632+1904 (L8) 


DENIS 0255-4700 (L9) 


SDSS 0423-0414 (TO) 


H 2 0-.ff 


0.74 ± 0.08 


0.71 


0.69 


0.64 


CH 4 -H 


0.98 ± 0.08 


1.08 


1.03 


0.99 


H 2 0-K 


0.76 ±0.14 


0.70 


0.69 


0.68 


ch 4 -a: 


0.89 ±0.14 


0.88 


0.84 


0.82 



Table 3. MS progenitor masses (Mms) and ages (tMs), WD cool- 
ing ages (tc), total age (t, ota i), and MS age as a fraction of WD 
cooling age for PHL 5038A. 







M WD (M Q ) 






0.57 


0.72 


0.87 


M MS (Me) 


2.1 


3.2 


4.4 


tMs (Gyr) 


1.3 


0.4 


0.2 


tc (Gyr) 


1.0 


1.5 


2.5 


Uotai (Gyr) 


2.3 


1.9 


2.7 


IMS / tc 


1.3 


0.3 


0.1 



aration from its primary of 55 AU. PHL 5038B has the potential 
of being used as a benchmark BD if its ages can be determined 
more accurately. A measurement of a trigonometric parallax for 
this system would enable us to more accurately constrain the ra- 
dius and, via models, the mass of the WD. This would then allow 
a better estimate for the total age of the system. 
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At an age of 1.9-2.7 Gyr (Table 3) the Lyon group models pre- 
dict a mass of around 60 Mj (Chabrier et al. 2000; Baraffe et 
al. 2001). Therefore, PHL5038B is almost certainly a bona fide 
substellar object. 

The mass luminosity relation for BDs is age dependent. 
Benchmark BDs can be used to test evolutionary models for sub- 
stellar objects if their ages can be satisfactorily constrained (e.g., 
Dupuy et al. 2009). One of the candidates for these studies is 
widely orbiting WD+BD binaries (Farihi, Becklin & Zuckerman 
2005; Pinfield et al. 2006) such as PHL 5038. The limiting factor 
in this system is its age calibration. Pinfield et al. (2006) require 
that a good age calibrator's progenitor lifetime be no more than 
10% of the WD cooling age (tMsAcX due to intrinsic uncertainty 
in calculating main sequence (MS) lifetimes. 

With r eff = 8000 + 100 K and logg = 8.2 + 0.1 the cool- 
ing age of PHL5038A is 1.0 - 2.5 Gyr (Fontaine, Brassard & 
Bergeron 2001). The mass of the progenitor to PHL5038A is 
estimated to be 3.2+1.1 M using the initial to final mass rela- 
tionship of Dobbie et al. (2006), which is valid for initial masses 
of > 1.6M (Kalirai et al. 2007). This indicates the WD de- 
scended from an A or late B-type star. An approximate main 
sequence lifetime can be calculated using the models of Girardi 
et al. (2000), giving 0.2 -1.3 Gyr. Therefore, the total age of the 
system is 1.9-2.7 Gyr. These ages were calculated for the lower 
and upper mass limits of PHL5038A and are listed in Table 3 
along with the ratio tMs/?c- This ratio ranges from 0.1 - 1.3 for 
a Mwd = 0.57 - 0.87 M and so only fulfills Pinfield's require- 
ment for a benchmark BD if the WD mass is as high as 0.87 M . 
At this point, we note that the masses for DA white dwarfs cooler 
than 12,000 K derived from spectroscopic fitting techniques may 
be over-estimated by up to 0.2 M , for reasons that are not fully 
understood (Kepler et al. 2007; Koester et al. 2008). 

6. Conclusions 

We have confirmed the presence of and spectroscopically typed 
a BD companion to the WD PHL5038A. PHL5038B was as- 
signed a spectral type of L8-L9 based on the measurement of 
standard spectral indices in the NIR, the first time this has been 
accomplished for a bona fide WD+BD binary. This spectral type 
is corroborated by the photometry of PHL5038B and by com- 
parison with the spectrum of an L8 field BD. The system is sit- 
uated 64 pc from Earth, and the BD has a projected orbital sep- 
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